Astronomy & Astrophysics manuscript no. just-sfr 


© ESO 2008 


February 1,2008 





The local star formation history of the thin disc derived from 

kinematic data 



o . 

O ■ Andreas Just 1 and Hartmut JahreiB 1 

a 
2 



Astronomisches Rechen-Institut at ZAH, MonchhofstraBe 12-14, 69120 Heidelberg, Germany 
Printed: February 1, 2008 

ABSTRACT 



.-?"*{ Context. We present an evolutionary disc model for the thin disc in the solar cylinder based on a continuous star formation history (SFR) and 
CLi a continuous dynamical heating (AVR) of the stellar subpopulations. 

^ Aims. We determine the star formation history of the thin disc in the solar vicinity. The self-consistent model of the vertical structure allows 
j.^ . predictions of the density, age, metallicity and velocity distribution of main sequence stars as a function of height above the midplane. 

Methods. The vertical distribution of the stellar subpopulations are calculated self-consistently in dynamical equilibrium. The SFR and AVR 
C^j , of the stellar subpopulations are determined by fitting the velocity distribution functions of main sequence stars. 

1 1 Results. We find a vertical disc model for the thin disc including the gas and dark matter component, which is consistent with the local 

j ' kinematics of main sequence stars and fulfils the known constraints on the surface densities and mass ratios. The SFR shows a maximum 
^ , lOGyr ago declining by a factor of 10 until present time. The velocity dispersion of the stellar subpopulations increase with age according 
^ ■ to a power law with index 0.375. Using the new scale heights leads to a best fit IMF with power-law indices of 1.5 below and 4.0 above 
I/-} | 1.6Mq, which has no kink around IMq. Including a thick disc component results in slight variations of the thin disc properties, but has a 
QQ ■ negligible influence on the SFR. A variety of predictions are made concerning the number density, age and metallicity distributions of stellar 
(f) subpopulations as a function of z above the galactic plane. 

Conclusions. The combination of kinematic data from Hipparcos and the finite lifetimes of main sequence stars allows strong constraints on 
^—-j . the structure and history of the local disc. A constant SFR can be ruled out. 

. Key words. Galaxy: solar neighbourhood - Galaxy: disk - Galaxy: structure - Galaxy: evolution - Galaxy: stellar content - Galaxy: kine- 
matics and dynamics 

• i-H , 

X: 

^ , 1. Introduction terminations of late type dwarfs. They apply stellar evolution 
- - and scale height corrections. The main result is the determina- 
The star formation history SFR(t) of the Milky Way disc is tion Qf enhanced star formation episodes over the lifetime of 
still not very well determined. The main reason for that is the ^ disc They exdude a const£mt SpR from chemical eyolu . 
lack of good age estimators with corresponding unbiased stel- ti(m modelf , Jn R & Q dg , a Fuente Marcofj ^ ^ sm 
lar samples. Especially samples selected by colour cuts or by a formation history for open star clusters were determined . They 
magnitude limit are biased with respect to the age distribution, found ^ ^ fiye episodefj of enhanced star formation in the 
because there is an age-metallicity relation due to the chemical ^ 2 Gyr since star duster contain Qnly a smaU percentage of 
enrichment of the disc. Therefore not even the famous Geneva- ^ disc ^ an extrapolation t0 the total (smooth ed) SFR is 
Copenhagen sample of 14.000 F and G stars (Nordstroem et not p 0SS j^i e 
al. |2004| l with well determined stellar properties including in- 
dividual age estimates can be used to derive the star formation In recent years the Hipparcos stars with precise parallaxes 
history directly by star counts. The most complete disc model and proper motions were used to determ ine the SFR with dif- 
based on star counts is that of Robin et al. (120031 1. They used ferent methods. Hernandez et al. ( 120001) determined the SFR 
a series of stellar subpopulations of different ages including over the last 3 Gyr using isochrone ages. They found a series of 
the AVR and the chemical enrichment but fixing the SFR to star formation episodes on top of an underlying smooth SFR. 
be constant. As a result, scale heights and surface densities of This result is complementary to our model, which gi ves the 
the stellar subpopulations are systematically smaller than in our slow evolution of the smoo thed SFR. In Binney et al. ( j2000| > 
model. Rocha-Pinto et al. (120001 1 used chromospheric age de- and in Cignoni et al. ( 120061 the scale height variation of main 

sequence stars were not taken into account. Therefore these 
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SFR. Binney et al. (2000) determined an age of the thin disc 
of 1 1 .2 Gyr (consistent with the 12 Gyr we are using) and both 
papers found an approximately constant age distribution, which 
is fully consistent with our model. 

We present the first disc model, which determines the 
SFR, the AVR and scale heights consistently. We use the 
Hipparcos stars and at the faint end the Catalogue of Nearby 
Stars (CNS4) to construct an evolutionary disc model includ- 
ing the kinematic information. We select the main sequence 
stars and divide these into a series of volume complete sub- 
samples from B to K type. In a kinematic model we derive the 
vertical velocity distribution functions /(| W\) of each subsam- 
ple, which depends mainly on the star formation history and 
the dynamical evolution of the disc described by the dynamical 
heating <Jw(t) (the age velocity dispersion relation AVR). The 
derived SFR depends on the normalized velocity distribution 
functions only. Therefore the model is essentially independent 
of the adopted IMF. In a final step we use the disc model to 
derive the local IM F from the volume complete subsamples 
applying consistently the lifetime and the scale height correc- 
tions. 



2. Disc Model 

In this section we describe the construction of the disc model 
and the determination of the properties of the stellar subsam- 
ples. For the self-consistent determination of the gravitational 
potential we take into account the stellar component, the gas 
component and the dark matter halo. 

We use a thin, self-gravitating disc composed of a sequence 
of stellar subpopulations according to the SFR and the AVR. 
Additionally a cold gas component and a Dark Matter halo are 
included for the gravitational forces. We take into account finite 
lifetimes of the stars and mass loss due to stellar evolution. The 
stellar remnants stay in their subpopulation, the expelled gas 
(stellar winds, PNs, and SNs) is mixed implicitly into the gas 
component. Since the lifetimes and mass loss depend on metal- 
licity, the metal enrichment with time is included. A standard 
IMF is used for the determination of the stellar mass fraction 
of the subpopulations with age. 

From large sets of template functions for the SFR and 
AVR we search for the best fitting pair of input functions. The 
total gravitational potential and the density profile of each age- 
bin are determined self-consistently assuming isothermal dis- 
tribution functions with velocity dispersion according to the 
AVR. The velocity distribution functions /(IW-'I) for main se- 
quence stars are calculated by a superposition of the Gaussians 
weighted by the local density contribution up to the lifetime of 
the stars. These distribution functions are simultaneously com- 
pared to the observed / (| W^Q in magnitude bins along the main 
sequence to judge the quality of the model. 

2.1. Self-gravitating disc 

The backbone of the disc is a self-gravitating vertical disc pro- 
file including the gas component in the thin disc approximation. 
In this approximation the Poisson-Equation is one-dimensional 



and in the case of a purely self-gravitating thin disc (i.e. no ex- 
ternal potential) the Poisson equation can be integrated leading 
to 

^ = 8ttG^ ■ (1) 

Therefore we model all gravitational components by a thin disc 
approximation. We include in the total potential $, the stellar 
component <3> s , the gas component $ g , and the dark matter halo 
contribution $h 

$(z) = $ s (z) + $ g (z) + $ h (z) . (2) 

In order to obtain the force of a spherical halo correctly in 
the thin disc approximation, we use a special approximation 
(see subsection |2.4| i. Since we construct the disc in dynamical 
equilibrium, the density of the sub-components are given as a 
function of the total potential p- : (<3>). The vertical distribution is 
given by the implicit function z(Q) via direct integration 
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The relative contribution of the stellar, the gaseous, and 
the DM-component to the surface density (up to \z\ = z max ) 
are given by the input parameters Q s ,Q g , Qh, which have to 
be iterated to fit all observational constraints. The disc model 
has two global scaling parameters to convert the normalized 
model to physical quantities. We fix the local stellar density to 
Pao = 0.039 Mq pc -3 , which is the best observed quantity for 
the local disc model (JahreiB et al. |19971 l. The second scaling 
parameter is the exponential scale height z s , which fixes via the 
shape correction the half-thickness of the stellar disc h e g, with 
the mass fractions Q s ,Q g , Qh the surface densities, and via the 
maximum velocity dispersion <j e the scaling of the velocity dis- 
tribution functions. 

2.2. Stellar disc 

The stellar component is composed by a sequence of isother- 
mal subpopulations characterized by the IMF, the chemical en- 
richment [Fe/H](t), the star formation history SFR(t), and 
the dynamical evolution described by the vertical velocity dis- 
persion aw(t)- Here r = t a — t is the time and t is the age 
of the subpopulation running back in time from the present 
time t a — 12 Gyr (which is the adopted age of the disc). 
We include mass loss due to stellar evolution and retain the 
stellar-dynamical mass fraction g(t) (stars + remnants) only. 
The mass lost by stellar winds, supernovae and planetary neb- 
ulae is mixed implicitly to the gas component. 

With the Jeans equation the vertical distribution of each 
isothermal subpopulation is given by 



p s ,i (z) = psoj exp 



-*(*) 



(4) 



where p s j is actually a 'density rate', the density per age bin, 
and er\v,j the velocity dispersion at age t y The connection to 
the SFR is given by the integral over z 
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with time Tj = t a — t; . The (half-)thickness h p (i; ) of the sub 
populations is defined by the midplane density p s0 j through 

gWSFIUn) 



PS ° J 2h p (t- 3 ) 
and can be calculated by 
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The total stellar density p s (z) and velocity dispersion <j s (z) are 
determined by 



Ps(z) 



Ps,j{z)dt 
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p s (z) then determines the potential & s (z) via the Poisson equa- 
tion. The total stellar surface density E s is connected to the 
integrated star formation So by the effective stellar-dynamical 
mass fraction g ff 



E s = J p s (z)dz = g e eS 



with 



So 



SFRdt and g eS 



J g(t)SFR{T)dt 



(10) 



(11) 



which includes luminous stars and stellar remnants. The local 
stellar density is given by 



PsO = 



gegSo 



(12) 



with thickness h e s f° r all stars. For a general shape of the den- 
sity profile h e s is different to z s , the exponential scale height 
well above the midplane. The effective exponential scale height 
z s of the stellar disc describing the exponential decline above 
the midplane is connected to the maximum velocity dispersion 
<7 e of the subpopulations and the total surface density Etot via 



C 7 Z„ — Cn 



27rGE t 



(13) 



where z c is the exponential scale height of an isothermal com- 
ponent above a disc with total surface density E to t- The shape 
correction factor C z is of order unity and is determined by the 
mean exponential scale height in the range z = (2 — 5) z e . 

The metallicity [Fc/H] affects the lifetimes, luminosities 
and colours of the stars and as a consequence also the mass loss 
of the subpopulations. In order to account for the influence of 
the metal enrichment we include a metal enrichment [Fc/H] (t) 
(see Fig. [TJ, which leads to a local metallicity distribution of 
late G dwarfs consistent with the observations (see Fig. [TTJ. 
The properties of the stars and the stellar subpopulations are 
determined by population synthesis calculations (see Sect. l2.5b - 

The properties of main sequence stars with lifetime r 
are determined by an appropriate weighted average over the 



age range. 
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the normalized density profile 
are given by 
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Density profiles and thicknesses of the final model are shown 
in Figs. [8] and [9] and the velocity dispersion of the main se- 
quence stars are compared in the lower panel of Fig. [5] with the 
observations. 

In the fitting procedure (Sect. 14. It a pair of star formation 
history and heating function is selected to derive the intrinsic 
structure of the disc. The star formation history and heating 
function of the final model are shown in Fig.Q] 

2.3. Gas component 

For the gas component we use a simple model to account for 
the gravitational potential of the gas. The vertical profile of the 
gas component, which is used for the gravitational force of the 
gas, is constructed dynamically like the stellar component. The 
gas distribution is modeled by distributing the gas with a con- 
stant rate over the velocity dispersion range o\y (t ) of the young 
stars up to a maximum age t g . By varying t g we force the scale 
height of the gas z g to the observed value of z g w 100 pc. The 
surface density of the gas is related to the stellar surface den- 
sity by the ratio Q g /Q s = E g /E s , which is chosen (together 
with the scale height h c g of the stars) to be consistent with the 
observed surface density of E g « 10 — 13 Mqpc~ 2 . 

2.4. Dark matter halo 

The halo does not fulfil the thin disc approximation. For a 
spherical halo we get the vertical component of the force to 
lowest order from 



d$ h _ GM R z 
~~dz~ ~ R? R 



(19) 



with r = R +z and the enclosed halo mass Mr inside radius 
R. Comparing this with the one-dimensional Poisson equation 
from the thin disc approximation (integrated over z near the 
midplane to lowest order for small z) 



— w AnGpoz 
dz 

we can use for the local halo density 

Mr 
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Fig. 1. The upper panel shows the SFR for the model and the 
infall rate of primordial gas derived from a local chemical evo- 
lution model (see Sect. 14.51 ). The thin line is the mean value 
of the SFR for comparison. The next panel gives the velocity 
dispersion of the stellar subpopulations The lower panel shows 
the adopted metal enrichment [Fe/H](t), which leads to a lo- 
cal metallicity distribution of late G dwarfs consistent with the 
observations (see Fig.Qj}. 



to be consistent with the spherical distribution. This value cor- 
responds exactly to the singular isothermal sphere. Therefore 
we can use the thin disc approximation also for the halo by us- 
ing the local halo density pho and the halo velocity dispersion 
(7h estimated from the rotation curve by adopting an isother- 
mal spherical halo. The effect of a cored halo, anisotropy and 
flattening, which would lead to small inconsistencies at large 
z, is neglected here. For other halo profiles correction factors 
would be necessary introducing some inconsistency in the halo 
profile description and leading to a different local halo density. 
The latter is be more important, because the effect of the halo 
potential on the disc is stronger than the adiabatic contraction 
of the halo in the disc potential. 

We use an isothermal halo component with Ch ~ 100 km/s. 
The halo surface density Eh is determined in the fitting proce- 
dure, because it influences not only the total surface density 
significantly, but also the shapes of the velocity distribution 
functions / (jM^I) of the stars. 



2.5. Stellar population synthesis 

For the determination of luminosities, main sequence lifetimes 
and mass loss we are not interpolating directly evolutionary 
tracks of a set of stellar masses and metallicities. In order 
to get a complete coverage of stellar masses we use instead 
the stellar population synthesis code PEGASE (Fioc & Rocca- 
Volmerange 1 19971 1 to calculate mass loss and luminosities of 
"pseudo" simple stellar populations (SSPs). This means that 
the PEGASE code is used to calculate the integrated luminosi- 
ties and colour indices for a stellar population created in a sin- 
gle star-burst at different time-steps. These SSPs are then used 
to assemble a stellar population with a given star formation his- 
tory, in the sense that the star formation history is assembled by 
a series of star-bursts. In this way we can assemble stellar pop- 
ulations for varying SFR and metal enrichment [Fc/H]. Our 
SSPs are modeled by a constant star formation rate with a du- 
ration of 25 Myr, the time resolution of the disc model. This is 
done for a set of different metallicities and intermediate values 
from the chemical enrichment are modeled by linear interpola- 
tion. 

The application of the PEGASE code is twofold. On one 
hand mass loss due to stellar winds, planetary nebulae and su- 
pernovae determines g(t), the mass fraction remaining in the 
stellar component as a function of age t. This depends on the 
IMF and the metal enrichment. We adopt a Scalo-like IMF 
(Scalo [T986i l given by 



dN oc M~ a dM 

1.25 0.08 < M/M© <1 

a = { 2.35 for 1 < M/M© < 2 
3.0 2 < M/M© < 100 



(22) 
(23) 



and five different metallicities ([Fe/H]= -1.23, -0.68, -0.38, 0.0, 
0.32), which are the input parameters into the code. Fig. |2] 
shows the mass loss for the different metallicities (thin full 
lines). The thin lines show the contribution from luminous mat- 
ter and stellar remnants for the set of input metallicities. The 
sum of both contributions for each metallicity vary by a few 
percent only and are therefore not shown. The thick full line 
shows g(t), the fraction of stellar mass in the present day stel- 
lar disc as a function of age including the chemical evolution. 
For the oldest age-bins about 40% of the stellar mass is lost by 
stellar evolution. The total fraction of stellar mass g c g = 0.644 
is indicated by the horizontal line. 

In the second application we determine the main sequence 
lifetimes and luminosities for stellar mass bins. Here we use the 
PEGASE code to compute V-band luminosities of isochrones 
in small mass bins. These are needed to estimate the mean 
main sequence lifetime as a function of the V-band luminos- 
ity My for the calculation of the velocity distribution functions 
/(|VF|) of these stars. We apply the PEGASE code to piece- 
wise constant IMFs for mass-bins with O.IMq and for a finer 
grid of metallicities ([Fe/H]= -0.8, -0.68, -0.5, -0.28, 0.0, 0.20, 
0.32). In Fig. [3] we show My for the different mass bins as 
a function of age. The lower panel gives the luminosity evolu- 
tion along the main sequence for the different metallicities (thin 
lines) and the age dependence of the stars in the disc model for 
M = 0.8 Mq demonstrating the significant variation over the 
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Fig. 2. Mass loss due to stellar evolution. The thin lines show 
the fraction of luminous mass and of remnants (upper and 
lower curve, respectively) for the set of metallicities used in 
the PEGASE code. The thick lines shows the total stellar mass 
fraction (luminous plus remnants) of the model taking into ac- 
count the age-metallicity relation. 

whole age range. The upper panels show the age dependent lu- 
minosities for all mass bins. These are again not the luminosity 
evolutions of the stars but the present day properties of the stel- 
lar population taking into account the age-metallicity relation. 
The vertical lines indicate the estimated mean main sequence 
lifetimes in the corresponding luminosity bins. The criteria to 
choose these lifetimes for the determination of the velocity dis- 
tribution functions /(|W^|) in the magnitude bins are discussed 
in Sect. 03] 

3. Observational Data 

Each stellar subpopulation, which is on average older than 
10 s yr, is in dynamical equilibrium with respect to the vertical 
distribution. Therefore they can be used as independent sam- 
ples to determine the total vertical gravitational potential <fi(z). 
The velocity distribution functions /(IM^I) of the samples de- 
pend on the age distribution of the stars. Therefore we use main 
sequence stars, where /(iM^j) is a function of the lifetime of 
the stars. In order to determine the vertical velocity distribution 
function / ( | W | ) in the solar neighbourhood we need kinemat- 
ically unbiased samples with space velocities. 

For the determination of the age velocity dispersion rela- 
tion we use the McCormick K and M dwarfs (Vvssotsky |19631 l. 
Detected by a spectroscopic survey they are free from kine- 
matic bias. Altogether 516 stars show reliable distances - al- 
most all from the Hipparcos Catalogue - and space velocity 
components. For a subsample of about 300 stars Wilson and 
Woolley (119701 1 estimated the Call emission intensity at the H 
and K lines in a relative scale allowing to construct six differ- 
ent age bins under the assumption of a constant star formation 
rate ( see JahreiB and Wielen, 1997). For each bin the vertical 
velocity dispersion is determined. The mean ages of the bins 
are rescaled from lOGyr to 12Gyr for the total disc age (see 
Fig. |5j. All these stars are also used to determine the velocity 
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Fig. 3. The lower panel shows the age dependence of the abso- 
lute luminosity My for stars with M = 0.8 Mq (averaged over 
0.1 Mq). Thin lines are the evolution for different metallic- 
ities [Fe/H] = -0.8; -0.68; -0.5; -0.28; 0.0; 0.20; 0.32 and 
the thick lines shows the present day luminosity of the stars as 
a function of age taking into account the age-metallicity rela- 
tion. 



distribution function /(|W^|) of stars with lifetime larger than 
12 Gyr (see Fig.[6]l. The number of stars in each individual age 
bin is too low to obtain reliable 

For the determination of /(iVFj) along the main sequence 
we use the Hipparcos stars with good space velocities supple- 
mented at the faint end down to My = 9.5 mag by stars from 
an updated version of the Catalogue of Nearby stars (CNS4, 
JahreiB et al. 1997), i.e. also most of the CNS4 data rely on 
Hipparcos results. For the determination of the space veloci- 
ties good distances, proper motions and radial velocities are 
required. This was achieved in combining the Hipparcos data 
with radial velocities originating from an unpublished compi- 
lation of the "best" radial velocities for the nearby stars, which 
was then extended to all Hipparcos stars. 

The selection process is a following. The visual abso- 
lute magnitude of the Hipparcos stars was determined from 
the visual magnitude and the trigonometric parallax given in 
the Hipparcos catalog. In the case of binaries resolved by 
Hipparcos the My of the brighter component was chosen cal- 
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Fig. 6. The normalized velocity distribution functions /(|W"|) for the different sub-samples given in Table Q] and for the 
McCormick stars (lower right). Symbols are the data and lines represent the model. 



culated from the combined magnitude and the magnitude dif- 
ference measured by Hipparcos. Then we select in the HR- 
diagram (Fig. [4J a rough regime along the main sequence in 
order to exclude most of the giants and white dwarfs. After 
determination of the best fit for the mean main sequence 
My, ms (B — V) (thick full line) all stars in the magnitude 
range My.ms ± 0.8 mag were selected. Only for the bright- 
est magnitude bin with My < 0.5 mag we use a colour cut 



B — V < 0.1 mag instead. Now the main sequence is divided 
into magnitude bins My < 0.5, 1±0.5, 6±0.5, 8±1.5mag. 

In order to avoid a kinematic bias we restrict the distances 
of the stars in each magnitude bin to be well within the com- 
pleteness limit of the Hipparcos catalogue determined by the 
magnitude limit V ~ 7.3 mag of the Hipparcos Survey. The 
properties of the resulting samples of main sequence stars 
are collected in Table Q] Column 1 lists the source catalogue. 
Column 2, 3 and 4 list the range in absolute visual magnitude, 
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Fig. 4. The HR-diagram shows all Hipparcos stars (grey dots) 
with ov/V < .10 The selected main sequence stars are over- 
plotted by larger black dots and the eight magnitude bins are 
indicated by the dotted lines. 



the selected distance limit, and the total number of stars avail- 
able, respectively. Column 5 and 6 list the number of stars re- 
moved due to poor parallaxes (ov/tt > .15) and unknown or 
poor radial velocities. Finally, in column 7 the remaining num- 
ber of main sequence stars with good space velocity compo- 
nents is given. 



Table 1. Complete samples of nearby main sequence stars 
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Fig. 5. Upper panel: Velocity dispersion as a function of age 
for the stellar subpopulations. Full symbols are the age groups 
of the McCormick K and M dwarfs. The asterisks reproduce 
the age bins of Fig. 31 in Nordstrom et al. (2004). The full 
line shows the heating function AVR of the final model. Lower 
panel: The symbols show the velocity dispersion of the stars 
along the main sequence (magnitude bins) as a function of 
mean lifetime. The full line gives the result of the final model. 



For the determination of the velocity distribution functions 
/(IV^I) we correct for the peculiar motion of the Sun using 
Wq = 7km/s. The resulting normalized distribution functions 
are shown in Fig. [6] with a binning of 5km/s in \ W\. 
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(*) n = number of stars with a-n/n > .15 



4. Properties of the disc 

The disc model has some free parameters additionally to the 
main input functions SFR and AVR. In all models we fix 
the local stellar density to p a o — 0.039Mqpc~ 3 from the 
CNS4 and the velocity dispersion of the Dark matter halo to 
Oh = 4<j e . For each pair of SFR and AVR we determine itera- 
tively the exponential scale height of the stellar disc z s , the rela- 
tive fractions Q s ,Qg, Qh of the surface density (at \z\ < z maK , 
where z max is indirectly determined by the depth of the poten- 
tial well), and the scale height of the gas component by choos- 
ing a maximum 'age' t g . The Q-values together with z s deter- 
mine the maximum velocity dispersion cr of the stars and the 
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surface densities of the components. A comparison with the 
observational constraints for the surface densities and veloc- 
ity dispersions restrict the exponetial scale hwight to z s =270pc 
with an uncertainty of 10%. In the final iteration the metal en- 
richment [Fe/H](t) is adapted to be consistent with the local 
metallicity distribution of G dwarfs. A list of model parameters 
and derived physical quantities is given in Table|2]together with 
the corresponding data of other authors. 

Table 2. Disc properties: mean and present day SFR (in 
Mo/pc 2 Gyr) are f° r the thin disc; local densities for stel- 
lar thin disc, gas, DM halo, thick disc (indices s, g, h, t; in 
Mq/pc 3 ); surface densities are below |z max | = 2.3 kpc or 
up to 0.35, 1.1 kpc (in M /pc 2 ); (half-)thickness of the thin 
disc h e s and exponential scale heights z are in pc; velocity dis- 
persions in km/s; The second column represents the thin disc 
model and the third column the values with the thick disc 1 
included. Columns 4 and 5 give some parameters from the lit- 
erature. 



quantity 


thin 
disc 


+ thick 
disc 


other 


sources 


(SFR) 


4.2 


3.9 






SFR P 


0.84 


0.79 






Po 


0.094 


0.092 


0.076 1 


0.098 2 


PsO 


0.039 


0.036 


0.045 1 


0.044 2 


PgO 


0.042 


0.041 


0.02 1 1 


0.050 2 


PhO 


0.013 


0.012 


0.01 1 




pto 




0.003 






s s 


32.6 


30.3 


34.4 2 




s s 


13.4 


13.3 


6 1 


13 2 






5.7 






^disc 


46.0 


49.3 


56 3 


48 4 


£(< 0.35) 


43.2 


42.4 


41 3 




E« 1.1) 


72.0 


71.5 


74 3 


71 4 




416 


418 






z B 


270 


270 






% 


100 


101 


140 1 




Zt 




745 






(To 


24.4 


24.3 


17.5 1 




Oh 


97.5 


97.1 


85 1 




Ot 




41.3 







1 Robin et al. |2003l 2 Holmberg & Flynn |2000l 3 Holmberg & Flynn 
[20041 4 Kuiken & Gilmore [T99T1 



In the next subsections we discuss the main input functions 
SFR and AVR and other properties of the disc model includ- 
ing a comparison with other work and predictions for future 
observations. 

4. 1. Star formation history and dynamical heating 

For a given stellar density profile there is a series of function 
pairs (SFR,AS/R) to construct a dynamical equilibrium model. 
Strong additional constraints originate in the correlation of the 
velocity distribution functions /(|VF|) of stars with different 
lifetimes. This leads to restrictions for the AVR and in turn to 



the SFR. The sequence of /(|W|) along the main sequence 
is a direct measure of the age distribution of stars in the local 
volume, which are converted to surface densities for connecting 
to the SFR. 

In the final model we use the star formation history (see 
Fig. ID 

SFR(t) = S ^^^ 2 with (24) 
a = 34.5 ; b = 4.2 ; t = 0.01 , (25) 

where the time t is in units of Gyr and So = 50.5 Mq/pc 2 is 
connected to the stellar surface density by Eq.[l0] The SFR de- 
creases by a factor of 10 from the maximum to the present day 
value, which corresponds to an e-folding timescale of 5 Gyr. 
The present day star formation rate with 20% of the average 
value of 4.1 Mq/pc 2 Gyr is relatively low. It is confirmed by 
the IMF considerations (see Sect. 14.71 ). 

For the dynamical heating function AVR (see Fig. [5) we 
use a power law 

o-w(f) = ( \ + f ,° ) with (26) 
\£a + to/ 

a = 0.375; t = 0.0767, (27) 

where t is the age in Gyr. The velocity dispersion of newly born 
stars is er p = 0.15cx o and tr = 24.4 km/s is the maximum ve- 
locity dispersion of the oldest stars with an age of t a = 12 Gyr. 
The power law index a — 0.375 is in the range of the classi- 
cal value of 1/2 from Wielen JT9771 ), of 0.47 (Nordstrom et al. 
12004b and the value of 1/3 derived by Binney et al. d2000b from 
proper motion measurements. 

4.2. Density profiles 

The vertical density profiles of the gas and of the stars are flat- 
tened at the galactic plane. A measure of the flattening is given 
by the ratio of the thickness and the exponential scale height. 

The corresponding stellar disc scale height is z s = 270 pc 
(from Eq. [13] with C z = 1.18), whereas the half-thickness h e g 
is 50% larger. 

For the gas, the thickness is 160pc compared to the scale 
height of lOOpc. The density profile and the surface density 
is consistent with the observed HI-profile (Dickey & Lockman 
( 1990) taking into account about 50% of H2 (Bronfman et al. 
(119881) within the large uncertainty range. 

Due to the gravitational potential of the disc, the local den- 
sity of the dark matter halo is 55% larger than the DM density 
at z = 2.3kpc and 18% larger than the mean DM density. 

The density profiles of main sequence stars differ in shape 
from the profiles of the subpopulations for a given age. The 
lower panel of Fig. [8] shows the normalized density profiles 
with lifetimes according to the samples used for the model. 
They are steeper than the corresponding profiles of the subpop- 
ulations with the same age and significantly shallower than the 
density profile using the mean age of the subpopulation. This 
difference is quantified by the thicknesses h p and h ms shown 
in the lower panel of Fig. [9] 
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Fig. 7. Vertical profiles of the disc. From top to bottom: Density 
profiles of the components; K z force, which is proportional to 
the cumulative surface density; velocity dispersion of the stellar 
component; V-band luminosity weighted mean metallicity of 
the stellar component. 
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0.1 
0.01 
0.001 
1e-04 
1e-05 




0.01 



0.001 



1e-04 



1e-05 



age =12Gyr 
=11Gyr 
=9.0Gyr 
=7.0Gyr 
=3.0Gyr 
=1.4Gyr 
0.8Gyr 
r 



0.5 



1.5 



z [kpc] 



lifetime >12Gyr 
=11Gyr 
=9.0Gyr 
=7.0Gyr 
=3.0Gyr 
=1 .4Gyr 
=0.8Gyr 
=0.3Gyr 




z [kpc] 

Fig. 8. Upper panel: Normalized density profiles of the stellar 
subpopulations as a function of age. Lower panel: Same for 
the main sequence stars as a function of lifetime. The set of 
lifetimes corresponds to the mean lifetimes of the subsamples 
used for the velocity distribution functions (see fig. [6). 



The density profiles derived by Holmberg and Flynn (2000) 
for two sub-samples of A stars (with < My < 1 mag) and 
early F stars (with 1 < My < 2.5 mag) are in full agree- 
ment with the corresponding profiles with lifetimes of 0.3 and 
0.8 Gyr in our model (see Fig. [8]). The K z force of their refer- 
ence model is also in good agreement with our results (see Fig. 

13. 

4.3. Age distributions 

The age distributions of the stellar sub-samples of MS stars 
depend on the lifetime and on the vertical distribution. Fig. [9] 
shows in the upper panel the age distribution of stars along the 
man sequence. Models of the star formation history, which cor- 
rect for the finite lifetime only and do not take into account the 
scale height variation, determine the local age distribution of K 
and M dwarfs. It varies by less than a factor of two around the 
mean value. Binney et al. (12000) and Cignoni et al. (2006) pro- 
pose a constant local age distribution in this sense. Taken the 
uncertainties of isochrone ages into account, this is consistent 
with our model. 



The effect of the scale height correction is demonstrated by 
including the SFR (thin line), which represents the age distri- 
bution of K and M dwarfs vertically integrated over the disc. 
Especially for young stars there is a significant overrepresen- 
tation in the solar neighbourhood. In order to quantify the cor- 
rection to the surface density, the lower panel of Fig. [9] shows 
the (half-)thickness h ms as a function of lifetime compared to 
the exponential scale height of the same populations z cxp and 
to the thickness of the stellar subpopulation as a function of age 
h p . 

The age distributions are a strong function of z above the 
plane. The middle panel of Fig.|9]shows the lack of young stars 
in steps of Az = 300 pc above the midplane. 

4.4. Kinematics of the disc 

The AVR can be observed only by kinematically unbiased stel- 
lar subsamples with direct age determinations. The upper panel 
of Fig. [5] shows the AVR of the model compared to two data 
sets. The circles are the McCormick K and M dwarfs with ages 
determined from H and K line strength fl 9971 and the asterisks 
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lifetime >12Gyr; z=0 

z=300pc 
z=600pc 
z=900pc 




age [Gyr] 



600 
500 



400 
300 




•a 200 

8 100 



6 4 
lifetime/age [Gyr] 

Fig. 9. The upper panel shows the normalized local age dis- 
tributions of main sequence stars with different lifetimes. The 
cumulative mass fraction of all stars and the SFR are also plot- 
ted. The middle panel shows the age distribution of K and M 
dwarfs and the cumulative stellar mass fraction with increas- 
ing height z above the plane. The lower panel shows the (half- 
thickness of the subpopulations as a function of age h p and of 
the main sequence stars as a function of lifetime h ms . For com- 
parison the exponential scaleheights z cxp for the main sequence 
stars are included. The horizontal line is the overall thickness 
of all stars h e g (Eq.fT2l) 



represent the F an G stars of Nordstroem et al. 2004 with good 
age determinations. 

The velocity dispersions of the subsamples of main se- 
quence stars with no age subdivision are determined by the 
weighted mean over the lifetime with the local age distribu- 
tion. The lower panel of Fig. [5] shows the excellent agreement 
of the model with the data from the nearby stars. Additionally 
the shape of the velocity distribution functions /(|W|) of the 
main sequence star samples are in good agreement with the 
model (Pig. [6j>. The shape of /(|W|) depends on the lifetime 
and also on the height z above the midplane. Fig. [10] shows 
the variation above the plane for K and M dwarfs (upper panel) 
and for F stars with My — 4 mag and a lifetime of 7 Gyr (lower 
panel). 



K, M dwarfs (M v >6); z=0 
z=300 pc 
z=600 pc 
z=900 pc 





20 30 

|W| [km/s] 

Fig. 10. The normalized velocity distribution functions /(| W|) 
change with increasing z (thick to thin lines) due to the vary- 
ing age mixture. The upper panel shows of K and M dwarfs 
with lifetimes larger than the age of the disc, whereas the lower 
panel is for late F dwarfs with lifetime 7 Gyr. 

4.5. Metallicity 

We constructed an analytic metal enrichment law [Fe/H](r), 
which reproduces the local metallicity distribution of late G 
dwarfs for the disc model with the final SFR and AVR. The lo- 
cal metallicity distribution is determined from the Copenhagen 
F and G star sample (Nordstrom et al. 2004) selecting all stars 
with masses 0.84 < M/Mq < 0.90 up to the completeness 
limit r < 40 pc. The lower mass limit of 0.84Mq was chosen 
in order to avoid an overrepresentation of metal poor stars. The 
metal enrichment law (Fig.Q]) is given by 



[Fe/H](r) 
[Fe/H] 



[Fe/H] +21g[l+pln(l + gT)] 
-0.7; p= 0.655; q = 0.833 



with 



(28) 
(29) 



with time r in Gyr. The initial metallicity is [Fe/H]o=-0.7 and 
the present day metallicity is [Fe/H] p =0.12. Before binning the 
theoretical distribution we add an observational scatter with 
FWHM=0. 165 dex. The comparison of the derived local metal- 
licity distribution for late G dwarfs with the data is shown in the 
upper panel of Fig. QT| The lower panels show the predicted 
metallicity distribution for late type stars (thick lines) and for 



A. Just, H. JahreiB: Local SFR of the thin disc 



11 



0.25 
0.2 

0.15 
0.1 

0.05 



N04 
model 



A 



A 



-1 -0.8 -0.6 -0.4 -0.2 0.2 

[Fe/H] 




-1 -0.8 -0.6 -0.4 -0.2 0.2 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




lifet>12.0Gyr 
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0.2 



Fig. 11. Upper panel: The metal enrichment law for the model 
(full line). The grey histogram shows the observations from the 
Copenhagen survey for late G dwarfs in the stellar mass range 
0.84 < M/Mq < 0.90. Lower panels: The metallicity distri- 
bution for early G dwarfs (lifetime 9 Gyr; thin lines) in compar- 
ison to the late G, K and M dwarfs (lifetime > 12 Gyr; thick 
lines) in the solar neighbourhood and 600 pc above the mid- 
plane. 



late F stars with lifetime 9 Gyr (thin lines) at the midplane and 
600 pc above the plane. 

The metal enrichment [Fe/H] and the SFR can be em- 
bedded in a local chemical evolution model. We proceed in 
the following way. We convert the metal abundance [Fe/H] 
to the oxygen abundance using the correlation of [O/Fe] and 
[Fe/H] from Reddy et al. (1993) by the linear approximation 
[O/H]=0. 375 [Fe/H]. For the oxygen enrichment we adopt in- 
stantaneous recycling and mixing to determine the infall rate 
of primordial gas. The mass loss from stellar evolution is 
taken into account. We start with a negligible initial amount of 
gas. Then the oxygen yield in solar units is given by y ox = 
(Z x) + -Zox.pHg/Es = 1.12 from the present day surface 
densities of gas and stars, the present day oxygen abundance 
\g(Z oxp ) = [0/H] p = 0.045 and the mean oxygen abundance 
(Z ox ) averaged over the star formation history. The infall rate 
and the surface density of gas and stars are shown in Fig. Q] 



0.1 



0.01 



0.001 



1e-04 



1e-05 



stars 
gas 

thick disc 1 :(41 , 6.7%, 1 .08, 745) 
2:(36, 7.8%, 1.41, 577) 
3:(29, 13.0%, 2.05, 380) 




z [kpc] 

Fig. 12. Self-consistent isothermal thick disc models with dif- 
ferent velocity dispersions. Stars, gas, thick disc 1 form a 
self-consistent model with parameters (cr t , Po,t/po,s+t, hi). 
Velocity dispersion and local mass fraction of the thick disc are 
input parameter and the power law index and exponential scale 
height are best fit parameter for a sech Q (z/ah t ) profile. For 
comparison two alternative cooler thick disc components are 
plotted. 



4.6. Thick disc 

Since the thick disc contributes a few percent to the local den- 
sity only, its properties are not well determined in the local disc 
model. Therefore we derive here as an example the density pro- 
file of an isothermal thick disc component and discuss the ef- 
fect on the thin disc. We split the local stellar density into a 
thin disc and a thick disc component. The local mass fraction 
of the thick disc is determined to equal the thin disc density at 
z « 1.2 kpc as indicated by K dwarf density profiles (Phleps et 
al. 12005). We choose the thick disc velocity dispersion and cal- 
culate the complete disc model including the new component. 
The parameters of all components have to be adjusted itera- 
tively, because the stellar disc profile (thin plus thick disc) has a 
different shape now. The resultant profiles for thin disc, gas and 
thick disc using er t = 41km/s and p Qit — 2.6 x 10 _3 Mq/pc 3 
corresponding to 6.7% of the local stellar density are shown in 
Fig.Q~2](thick lines). The density profile of the thick disc can be 
fitted by a sech Q (z/a/it) profile to better than 3%. The correc- 
tions of the thin disc parameters due to the thick disc is of the 
order of 1%. Only the surface density of thin and thick disc is 
10% larger with a corresponding reduction of the Dark matter 
(see Table|2]i. For two alternative models with smaller velocity 
dispersion of the thick disc the density profiles can also be fit- 
ted very well by the sech Q law. The parameters are also given 
inFig.E] 

4. 7. initial mass function 

For the determination of the IMF from the local luminosity 
function, the luminosity function is first converted to a mass 
function dN — f(M)dM using the transformation formulae 
of Henry & Mc Carthy |l993t l, corrected in Henry et al. d 19991 , 
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and extended to bright stars according to Schmidt-Kaler in 
Landolt-Bornstein (1982). The result is shown in the lower 
panel of Fig. [T3] and compared to the PDMF given in Kroupa 
et al. (|19931 l. The star numbers are normalized to a sphere with 
radius R=20pc. In order to determine the correction factor con- 
verting the IMF to the local PDMF due to the finite lifetime 
and the vertical thickness h ms , the main sequence lifetime is 
needed. We use the lifetimes estimated from the isochrone evo- 
lution shown in Fig. [3] for the different magnitude bins, where 
the metal enrichment and the relative weighting due to the in- 
creasing number of stars with decreasing mass in the mass in- 
terval is taken into account. These data are shown in the upper 
panel of Fig. Qj] with a comparison of the lifetimes determined 
directly from evolutionary tracks of Girardi et al. (2004) and 
with the analytic fitting formula of Eggleton et al. ( 1989). Any 
systematic variation of the lifetimes result in significant correc- 
tions to the correction factor and therefore the IMF. This is the 
most uncertain input to the IMF determination. The correction 
factors in the solar neighbourhood (z=0 thick line) and for pop- 
ulations above the plane are shown in the middle panel of Fig. 

El 

The used Scalo IMF (Eq. l22b is consistent with the ob- 
served PDMF. But the kink near 1 Mq seems artificial, because 
it is near that mass, where the correction factor due to the finite 
lifetime starts to apply. Therefore we determined a new IMF by 
fitting power laws in two mass regimes only. The best fit is 

dN oc M~ a &M with (30) 

f 1.49 ± 0.09 0.08 < M/M < 1.61 
a = < for 

{ 3.95 ± 0.33 1.61 < M/M Q < 6.0 

No = 359/M at M/M© = 1.61 , 

where Nq gives the normalization in the 20 pc sphere. The 
back-reaction of the corrected IMF to the disc model via the 
mass loss is very small and not included in the model. 

The strongest constraints on the IMF at high masses and 
the present day SFR is the observed number of A stars in 
the solar neighbourhood. Since the bright stars with My < 
0.5 mag are observed in a sphere with a radius of 200 pc, the 
sample size is a direct measure of the local surface density. 
Therefore the conversion from the IMF to the mean SFR in 
the last few 100 Myr depends only on the lifetime of the stars. 
That means, a higher present day SFR requires a shorter life- 
time for the A stars or an even steeper IMF. 

The position of the Sun is probably by about 20 pc above 
the midplane (Humphreys & Larsen[l995). The local density 
of the subpopulations with small scale height are significantly 
smaller than the midplane density (My < 1.5 mag). We cor- 
rected for that implicitly, since we determined the midplane 
density of these magnitude bins for the PDMF using spheres 
with large radii (see TableQ}, where the offset can be neglected. 

5. Summary 

We presented an evolutionary disc model for the thin disc in 
the solar cylinder based on a continuous star formation history 
(SFR) and a continuous dynamical heating (AVR) of the stellar 




-3 1 1 1 1 1 1 

-1 -0.5 0.5 1 

lg(M/M sun ) 

Fig. 13. Upper panel: Stellar mass as function of lifetime for 
the present day population. Symbols are the values used (from 
Fig. 0), the lines are estimated values from evolution tracks. 
Middle panel: Correction factor (=IMF/PDMF), which must be 
applied to the IMF to derive the local star counts of stars with 
corresponding lifetime relative to low mass stars with lifetime 
exceeding the disc age. Lower panel: Open circles show the 
PDMF and the IMF of our sample from CNS4 and Hipparcos. 
Crosses are the PDMF of Kroupa et al. ( 119931 ) and the IMF by 
correcting with the same dilution factor. The dashed lines are 
the Scalo IMF and PDMF used in the model and the full lines 
give the new best fit for the derived IMF (open circles). 

subpopulations. The vertical distribution of the stellar subpop- 
ulations are calculated self-consistently in dynamical equilib- 
rium. The SFR and AVR of the stellar subpopulations are de- 
termined by fitting the velocity distribution functions of main 
sequence stars. A chemical evolution model with reasonable 
gas infall rate is included, which reproduces the local [Fe/H] 
distribution of G dwarfs. 

We found a vertical disc model for the thin disc including 
the gas and dark matter component, which is consistent with the 
local kinematics of main sequence stars and fulfils the known 
constraints on the surface densities and mass ratios. The SFR 
shows a maximum lOGyr ago declining by a factor of 10 until 
present time corresponding to an e-folding timescale of 5 Gyr. 
The velocity dispersion on the upper main sequence depends 
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on the lifetime of the stars, which is derived from the AVR. For 
the AVR we find a power law with index of 0.375. Applying 
the stellar lifetime and the new scale height corrections to the 
PDMF results in a best fit IMF with power-law indices of 1.5 
below and 4.0 above 1.6 Mq, which has no kink around 1 Mq. 

Including a thick disc component consistently lead to slight 
variations of the thin disc properties, but a negligible influence 
on the SFR. A variety of predictions were made concerning 
the number density, age and metallicity distributions of stellar 
subpopulations as a function of z above the galactic plane. 
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